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Abstract 

We study low-lying states of the XY and Heisenberg antiferromagnets on a triangular 
" lattice to clarify whether spontaneous symmetry breaking occurs at T = in the ther- 

modynamic limit. Approximate forms of low-lying states are proposed, in which degrees 
I , of freedom of the sublattice magnetization and of the chirality are separated. It is shown 

,-0 , that low-lying states can be accurately described with the present approximation. It was 

D ■ 

["t i . argued that low-lying states play an important role in symmetry breaking. With help of 

this approximation, we discuss the contribution of low-lying states to symmetry break- 
ing of two types, namely creation of the spontaneous sublattice magnetization and the 
spontaneous chirality. Furthermore, to show an evidence for the occurrence of the sym- 
. metry breaking, we numerically study the low-lying states of finite systems of the XY and 

Heisenberg antiferromagnets. It is found that the necessary conditions for the symmetry 

CNJ ■ breaking to occur are satisfied in these models. 
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o ! 1 Introduction 

o ■ 

Many authors have studied ground-state properties of quantum antiferromagnets, e.g., the 
^ ! Heisenberg and XXZ models, on a triangular lattice. The Hamiltonian of the XXZ antiferro- 
magnet is given by 

h = jJ2 ( s i s j + s i s ] + ASfs;) (i.i) 

where the summation runs over all the nearest-neighbour pairs. Most of the studies are based 
on one of the following two viewpoints. The first one is that the system is disordered. Anderson 
proposed the RVB stated and Kalmeyer and Laughlin discussed a spin-liquid state® from 
this aspect. The second one is that the system has a long-range order with the so-called 
120° structure, though quantum fluctuations reduce the sublattice magnetization. Trial wave 
functions which have the 120° structure are discussed by Miyashita®, and by Huse and Elser^. 
The spin- wave expansions® BUS and the series expansion© suggested the existence of a long- 
range order. The present author and Suzuki® presented the exact ground states of the XXZ 
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model with A = —0.5, most of which have the 120° structures. We also studied the XXZ 
antiferromagnet in the region —0.5 < A < 1 using the spin-wave expansion and showed that 
the quantum fluctuations are enhanced, as A becomes large. The quantum effect is strong in 
the Heisenberg model. The estimates of the ground-state energy in the studies based on these 
pictures are almost equal to each other in the Heisenberg model. Thus it is difficult to judge 
which picture is correct. 

Apart from the above arguments, several authors© 0' & studied ground states in finite 
systems using the exact-diagonalization method. There are two papers© i-* which reported 
the calculations up to N = 36. Conclusions do not coincide with each other. Bernu, Lhuillier 
and Pierre® reported that the sublattice order is almost equal to that obtained by the spin- 
wave expansion. Leung and Runge's result® is, however, contradictory to the existence of the 
sublattice order. These two groups fit the data in different ways. 

Understanding of symmetry breaking has progressed by studies of low-lying excited states. 
In quantum antiferromagnets the ground state of a finite system is, in many cases, unique and 
symmetric. In the thermodynamic limit, symmetry is broken by making a linear combination of 
the ground state and low-lying excited states. Horsch and von der Linden© first showed that, 
if the Heisenberg antiferromagnet has a Neel order, there exists a low-lying excited state whose 
energy converges to that of the ground state. Koma and Tasaki© rigorously proved that, if 
the system has a long-range order, low-lying states of growing numbers converge to the ground 
state in the thermodynamic limit. Several numerical calculations showed the existence of these 
low- lying excited states in the spin 1/2 Heisenberg antiferromagnet on a square lattice© and 
on a triangular lattice©. 

Bernu et a/.© determined the quantum numbers of low-lying states which construct the 
ordered infinite-volume ground states in the Heisenberg antiferromagnet on the triangular lat- 
tice. They studied finite systems and found the existence of a whole set of low-lying states. 
Azaria et al. © studied the finite-size dependence of these low-lying states. They found that 
these states satisfy the scaling property which indicates the existence of a long-range order. 

In this paper we consider the XY and Heisenberg antiferromagnets on the triangular lat- 
tice. We study low-lying states to clarify whether spontaneous symmetry-breaking occurs at 
T = in the thermodynamic limit. We give approximate forms of low-lying states, in which 
the degrees of freedom of the sublattice magnetization and the chirality are separated. We 
show that the low-lying states can be accurately described with the present approximation. 
With help of our approximation we discuss how rearrangements of the low-lying states bring on 
symmetry breaking. In quantum antiferromagnets on the triangular lattice, there are two types 
of symmetry breaking; In the Heisenberg antiferromagnet both the sublattice magnetization 
and the chirality relate to the breaking of the 0(3) symmetry and in the XY antiferromagnet 
the 0(2)xZ2(chiral) symmetry can be broken. These symmetry breakings relate to each other. 
Using our approximation we can understand breakdown of each symmetry independently. Fur- 
thermore, studying the low-lying states of finite systems with the numerical diagonalization 
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method, we obtain an evidence of occurrence of the symmetry breaking. In the XY model 
the necessary conditions for the symmetry to break are fully satisfied. We conclude that there 
exists long-range order in the thermodynamic limit. In the Heisenberg model, since the system 
size studied here is too small, we can not still obtain a definite conclusion about the existence of 
long-range order, but there are many evidences which support the existence of the long-range 



order, as discussed in refs. 12 and 16 



The relation between the long-range order parameter a and the spontaneous sublattice 
magnetization m is also argued, though this is not the main part of the present paper. It has 
been discussed© that, for bipartite systems, the relation is given by m = y/3a in the Heisenberg 
antiferromagnet and by m = \[2a in the XY antiferromagnet. In the antiferromagnets on the 
triangular lattice, however, the factor of the relation becomes y/2 times as large as the above. 
We obtain m = ^/6a for the Heisenberg model and m = 2a for the XY model. 

Contents of this paper are the following. Approximate low-lying states are formally given 
and properties of these states are discussed in Section 0. Some results for the low-lying states 
of ferromagnets, which are used in Section ||] to discuss the properties of the approximation, 
are shown in Appendix A. The mechanism of symmetry breaking is discussed in Section |3]. In 
Section |4] we show an evidence of the symmetry breaking studying finite systems numerically. 
Section |5| contains summary and discussions. The relations between the symmetry breaking 
and the long-range order parameter are discussed in Appendix B. 



2 An approximation for low-lying states 

We study properties of the low-lying states of the quantum antiferromagnets on a finite- volume 
triangular lattice, using approximate states. In antiferromagnets on finite systems the ground 
state is symmetric. Symmetry breaking in an infinite-volume limit is discussed in the next 
section. 

We consider the quantum XXZ antiferromagnets on the finite- volume lattice A with a peri- 
odic boundary condition. The size of the system is N, where N = 3m with an integer m. The 
Hamiltonian is given by 

H = J £ (SfSf + SfSf + AS" S') , (2.1) 

(i,j}eA 

where the summation runs over all the nearest-neighbor sites and the symbol A denotes the 
anisotropy of the ^-component interactions. The models in the region —0.5 < A < 1 are 
considered. 

We concentrate on low-lying states which have sublattice-translational invariance and C%y- 
invariance. (We denote the point group of the 120° rotation and reflection of the lattice as C%y.) 
Bernu et al. reported that there are many low-lying states of this type. Here we classify 
the states according to eigenvalues of the 60° rotation (C 6 ). We call the class of C 6 -symmetric 
states as the type a and the class of C 6 -antisymmetric states as the type (3. 



3 



As an approximation for the lowest states of the type a and of the type (3 in the S^ otal = n 
(or n+ 1/2) subspace, we consider the following states: 

lna) \\(U + W)\Fn)\y \\(U-W)\Fn)\\ ^ 

for n = 0, ±1, ±2, . . ., where the symbol U denotes the unitary operator 

C/ = ex P (z^E^-4E^) (2-3) 

6 ieB 6 iec 

and \\(U ± 17+) |F?z> || = (Fn\(U ± U ] f\Fn) 1 / 2 . The state |Fn) denotes the lowest state in the 
S tot&\ = n (or n + 1/2) subspace of the "ferromagnetic" XXZ model on A which is defined by 
the Hamiltonian 

H F = -J {S?S* + SfSy -2AS*S*). (2.4) 

(ij)6A 

Here the parameter A is set equal to that in (|2.1| ). We denote the eigenvalue of Hp for the 
state \Fn) as Ef n . 

We first discuss the properties of the above approximate states ( |2.2|) . After that we verify 
that this approximation is good. 

Here we discuss the meaning of the referred model ( |2.4j ). The states \na) and \nf3) are linear 
combinations of U\Fn) and U'\Fn). The unitary operator U transforms the Hamiltonian ( |2.1|) 
into the form 

UHtf = — Jj2(S?Sj + S'fS] - 2AS?S*) - —J ( s ? s j - s t s j)- ( 2 - 5 ) 

Here the symbol % — > j goes from the sublattice A to B, B to C, and C to A. The first term in 
( |2.5|) gives the Hamiltonian (|2.4|) . Thus the present approximation corresponds to neglecting 
the second term in (|2.5| ). 

In the thermodynamic limit of the ferromagnet ( |2.4p , only 0(2) symmetry breaking can 
take place, since it has no frustration. The classical limit of the ferromagnet (|2.4|) in the region 
—0.5 < A < 1 has ferromagnetically ordered ground states, in which all spins are lying in 
the XY plane®. The spin 1/2 ferromagnetic XY model (A = 0) on the triangular lattice 
has been studied by using the exact-diagonalization method© and the spin-wave theory© §. 
The results indicate the existence of the ferromagnetic long-range order. As A increases in the 
ferromagnet (|2.4|), quantum fluctuations are enhanced and the magnetization is reduced. 

As shown in Appendix A, several exact results can be obtained for the ferromagnet ( |2.4| ) 
on the finite-volume lattice A. The ground state is unique and it exists in the S^ otal = 
subspace. Then it has 0(2)-rotational invariance. (When N is an odd number, there exists 
trivial degeneracy; the states |F0) and |F— 1) are degenerate.) The lowest state in the S* otal = n 
subspace, which we denote as |Fn), is unique. The coefficients of \Fn) in the basis which are 
eigenstates of { Sf } are nonnegative. All states { \Fn) } are translationally invariant and 
C 6 y-invariant. (See Appendix A.) 
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Table 1: Overlaps between the ground state |0) and the approximate state |0a), and between 
the first excited state |1) and the approximate one |0/3). 



oize jv 


Q 19 91 

o y rz zi 


XY K°l°«)l 2 
\{iW)\ 2 


1.0 0.9797 0.8472 0.9181 
1.0 0.9797 0.9988 0.9181 


|(0|0a)| 2 

Heisenberg 

|(1|0/3>| 2 


1.0 0.7881 0.7997 0.4643 
1.0 0.7881 0.9871 0.4643 



Using these rigorous results for the low- lying states of the ferromagnet, we can determine 
the spatial quantum numbers of the approximate states. We have verified in the small clusters 



that the spatial quantum numbers of the approximate states ( |2.2| ) are the same as those of 
the true states. Bernu et a/.© investigated the spatial quantum numbers of the low- lying 
states for the Heisenberg antiferromagnet and classified the low-lying states into three classes; 
Ti (states with k = and even under inversion), T 2 (states with k = and odd under inversion) 
and Instates with k = ±(4vr/3,0)). When 5 t z otal (= n) = 31 or S t 2 otal (= n + 1/2) = 3/ + 3/2 
with an integer Z, \na) and \n/3) belong to Ti and I^, respectively. Otherwise, both \na) and 
\n(3) belong to IV As we show in the next section, our classification is useful in discussing the 
symmetry breaking about the sublattice magnetization and the chirality separately. 

We now discuss accuracy of our approximation. We compare the approximate low-lying 
states with the true ones, numerically diagonalizing finite systems. First we note that, as we 
have already mentioned above, the spatial quantum numbers of our approximate low-lying 
states are equal to those of the true states: They are invariant under any sublattice translation 
and any translation of C^y, and they have the same momentum as the true states. In the true 
system there exist low-lying states of two types corresponding to the types a and (3. Secondly 
the approximate low-lying states have accurate expectation values of the energy, which are at 
most 2% higher than the exact values in the XY model and at most 7% higher in the Heisenberg 
model. Finally we calculated overlaps between the exact low-lying states and the approximate 
ones. The results are shown in Table [1]. In the XY model the approximate states have more 
than 90% overlaps with the exact states and in the Heisenberg model more than 70%. These 
situations are the same for the low-lying states in the S* otgl = 1 and S% otal = 2 subspaces. Thus 
we find that the exact low-lying spectrum has the two- fold structures, as shown in Fig. 1, and 
that the true low-lying states can be accurately described with the present approximation. 

The approximate ground state, |0a), is a generalization of the trial state of Betts and 
Miyashita® for the XY antiferromagnet. In the case of A = —0.5, |0a) and 1 0/3) belong to 
the exact ground states, which were presented by the author and Suzuki®. This fact partially 
gives the reason for accuracy of the state U\F0). 
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An advantage of our approximation is that the degrees of freedom of the ^(chiral) symmetry 
and the 0(2) symmetry are separated in it. Then we can discuss breakdown of each symmetry 
independently. The parts of the unitary operator, namely (U ± W), can describe the chiral 
symmetry breaking and the parts of the low- lying states of the ferromagnet, namely |Fn) 
(n = 0, ±1, ±2, . . .), can display the breakdown of the 0(2) symmetry. 

Using this approximation, we can estimate the expectation values of the energy. The ground- 
state energy is calculated as 

/n , WIn v £fo + 4ff f0 (F0|[/|F0) - 6A J(F0\U(E {hj) gg)|FQ) 

(Oa\H Oct) = 7— lrT .^ . — , 2.6) 

x 1 1 ' 2 + 2(F0|£/|F0) v ; 

where we have used the relations UHW + WHU = H F and UHU = H F W + WH F - 
3AJUiE(ij)SfSj. The terms (F0|Z7|F0) and (F0|£/(E(i,j> S?S?)\F0) correspond to the tran- 
sition probabilities from the ferromagnetic ground state, |F0), to the antiferromagnetic state, 
U\F0). In the thermodynamic limit, these values are vanishing. Thus the ground-state energy 
is estimated as 

{0a\H\0a) ~ (2.7) 

Next we discuss the energy gaps. The energy gap between the ground state (the type a) and 
the lowest state of the type /3 is almost vanishing as 

y , -3£fo(F0|[/|F0) + 6AJ(F0\U(E(ii)S"S")\F0) 
(OWm - (0a\H\0a) = \ _ p^p^ ,J> J ~ 

~ 0. (2.8) 

The energy gap between the ground state of S£ ot3l = and the lowest state of S^ otal = 1 is 
estimated as 

(la\H\la) - (0a\H\0a) ~ -{E n - E f0 ). (2.9) 

Relations between these energy gaps and the symmetry breaking of two types, namely creation 
of the spontaneous chirality and of the spontaneous sublattice magnetization, are discussed in 
Section || 

We give a remark for the case of the Heisenberg antiferro magnet. Since the model is 
isotropic, a slight modification is necessary for the approximate ground state, |0a(A=i))- We 
only give the form of an approximate ground state as follows: 

\0a') = —J dnR(n)\0a (A=1) ) /\\ — J dQR{n)\0a {A=1) )\\, (2.10) 

where Q is the unit vector with the spherical coordinates (8, if) and 

J2(0) =exp(ip£S?)«p(»*ESj)- (2.11) 

i j 

We used the state |0a(A=i)) to compare our approximation with the true ground state of the 
Heisenberg model, for it is not easy to use |0a') in numerical calculations. 
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3 Contribution of the low-lying states to the symmetry 
breaking 

Several authors© @0§) have discussed the mechanism of symmetry breaking in quantum 
Heisenberg antiferromagnets. It has been argued that, as the system size is enlarged, many 
low-lying states converge to the ground states and that linear combinations of these states 
make the symmetry breaking occur. In this section we study the roles played by the low-lying 
states of the two types in the symmetry breaking. Here we assume that the ground states of 
an infinite system have a long-range order and that the symmetry breaks down. In Section |] 
we show an evidence of occurrence of the symmetry breaking. 

In antiferromagnets on a triangular lattice, there are two types of symmetry breaking. Then 
the mechanism of symmetry breaking is complex. In the XY model, or in the XXZ model for 
—0.5 < A < 1, the 0(2) symmetry and the ^(chiral) symmetry can break independently. 
In the Heisenberg model, both the spontaneous sublattice magnetization and the spontaneous 
chirality correspond to the breakdown of the 0(3) symmetry. 

As shown in Section |2|, the low-lying states on a finite-volume triangular lattice have the 
two-fold structures. The purpose of this section is to clarify the relations between the symmetry 
breaking of the two types and the low-lying states of the types a and (3. Our approximation 
(|2.2|) helps us to understand the symmetry breaking of the two types independently, since the 
degrees of freedom of the sublattice magnetization and of the chirality are separated in it. 

In an infinite system we can define four types of ground states, which are classified by 
existence of the spontaneous sublattice magnetization or the spontaneous chiral order. Here 
we show the definitions of the four types. 

(1) Symmetric states (Mixed states), 

<• ■ •>! = Km Km ^rr^^ = • •)■ (3-1) 
x 71 7VToc/3Too Tr[exp(-/3#)] v ' y ' 

(2) States in which only the spontaneous chirality exists, 

x 11 bio jvtoo /3Too Tr[exp(-/3(# - BQ Z ))} v ; 

where Q z denotes the z-component of the chirality order-parameter operator 

Q4E(^ S^. (3.3) 

^ 3 <i-i> 



The summation is defined in the same way as in eq. 
(3) States in which only the spontaneous magnetization exists, 

/ v =r V r Tr[- ■ ■ exp{-/3(H - B(U SfW + W SfU))}] 

{ )3 b^nZ™ Tr[exp{-P(H - B(UEi.SfW + WEiSfU))}} ' l ' J 
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(4) Pure states, 

( ) - lim lim lim ^ - - ^~^ H Z gg? ^ S ? U ^ - Q( ) ( 35 ) 

Here the equilibrium states are defined as functionals of expectation values on an operator 
space. (This is familiar in studies of infinite- volume systems.) In the following, we show, using 
the approximate states (|2.2|), how each ground state is constructed of low- lying states. 

(1) Symmetric states (Mixed states). By taking the thermodynamic limit of the finite- volume 
ground state, we get the symmetric ground state ( |3.1|) . No symmetry is broken in it. 

(2) States in which only the spontaneous chirality exists. We obtain the state (• • by taking 
the thermodynamic limit under an infinitesimal effective field which is conjugate to the chiral 
order parameter. 

An approximation of this state is given by |2) = U\F0). There remains the 0(2) invariance 
in it. Both Miyashita's trial function© and the variational function by Huse and Elser§ belong 
to this type of ground states. 

In the XF-like model, the state C/|F0) is constructed by making a linear combination of 
\0a) and |0/3), 

{|0a>, |0/3)} >{U\F0), U^\F0)}. (3.6) 

Thus with help of our approximation we find that, if the ground state (the type a) and the lowest 
state of the type /3 become degenerate in the infinite-volume limit and if they are rearranged 
between themselves, the Z 2 symmetry of the chirality is broken. 

In the Heisenberg model, the degree of freedom of the chirality is continuous and the ground 
state is isotropic, as mentioned in fl2.10| ). Then two states are not enough to break entirely 



the symmetry of the chirality. Koma and Tasaki© proposed an approximate state which 
describes breakdown of continuous symmetry and they argued that the thermodynamic limit 
of this approximate state is a pure state, in which continuous symmetry is broken. By applying 
their argument to the present case, an approximate state which has fully-ordered spontaneous 
chirality can be constructed in the form 

(3.7) 




where Q = Q x ±iQ y and |$gs) denotes the ground state. The number k is of o(N). This state 
has a nonvanishing expectation value of the operator Q x . By rotating all the spins through 
the angle 7r/2 about the y axis, we obtain an approximation for the state (• ■ -)2- Following the 
arguments by Koma and Tasaki©, the infinite-volume limit of this approximate state become 
(• • -) 2 . The states (Q ± ) 2m |$ G s) belong to the type a and (Q ± ) 2m+1 |$ G s) to (3. Thus the state 
(• • -) 2 is constructed by a linear combination of low-lying states of the type a and of the type 
(3- 
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In our approximation the expectation value of the spontaneous chirality is calculated as 



(2|Q*|2> = ^ + |^F0|f/|F0>-2A(F0|£^|F0> 

(id) 



L[[) 2A(F0| SfS*\F0) (3.8) 



J 



(id) 



1 3 

where the relations UHU^ = -^H-p — -JQ and H = —H F + 3JAJ2(ij) S*Sj are used and E i0 
denotes the ground state energy of the ferromagnet. As A is increased from —0.5, the first 
term in ( |3.8|) becomes large by quantum effects. The second term in ( |3.8|) , (F0| J2{i,j) S*Sj\F0), 
behaves as follows: At A = —0.5 this term is vanishing, since there is no correlation between up 
and down spins in the ground state |F0). As the parameter A increases, nearest-neighbor pairs 
of up and down spins are favored. Then the term becomes large. Thus at A = the chirality 
is enhanced by the quantum effect of the first term and near A = 1 the value is reduced by 
contribution of the second term. 

(3) States in which only the spontaneous magnetization exists. In the state (• • -) 3 there exists 
the spontaneous magnetization, though there exists no spontaneous chirality. 

In the XF-like model, the 0(2)-symmetry breaking occurs in the following way. First 
we consider the case of the ferromagnetic XX Z model ( |2.4|) . As shown in Appendix A, the 
finite- volume ground state of the ferromagnet is unique and it has 0(2)-rotational invariance. 
It was shown in ref. 14 that the low-lying states \Fn) (n = ±1, ±2, . . . , ±o(N)) converge 
to the ground state and pure infinite-volume ground states are constructed by taking linear 
combinations of the low-lying states. We denote one of the ferromagnetically ordered ground 
states as |F-order), where (F-orderlS"? |F-order) ^ and (F-order| S'f |F-order) = 0. In the 
case of the antiferromagnet, with help of the approximate low-lying states, we can understand 
the mechanism of the 0(2)-symmetry breaking on the analogy of the ferromagnet. Under the 
effective field B in (|3.4| ), the low- lying states of the type a have an energy lower than those of 
f3 have. An approximate state for (■ • -)s is given by (U + £/')|F-order), which is constructed by 
making a linear combination of (U + W)\Fn) (n = 0, ±1, ±2, . . .). In the same way, by using 
low-lying states of the type (3, another state in which only the spontaneous magnetization exists 
is constructed as (U — L^)|F-order). Thus we find that, if the low-lying states of the same type 
become degenerate in the infinite-volume limit and if rearrangements occur between them, the 
breakdown of the 0(2) symmetry occurs. 

For the Heisenberg model, by applying the approximation of Koma and Tasaki© to this 
model, an approximate state of (■ ■ is written as 

(3.9) 




y/2k + 1 

where 0+ = UEiSfU* + U ] Y. % SfU. All the low-lying states, (0+) m |$ GS ) and (0-) m |$ G s), 
belong to the type a. Thus We can obtain the state (• • by taking a linear-combination of 
the low-lying states of the same type. 
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(4) Pure states. The pure infinite-volume ground states have the same structure as the 
ground states of the classical model. The 0(2) symmetry and Zi symmetry are broken in the 
XY model, and the 0(3) symmetry of the sublattice magnetization and that of the chirality 
are broken in the Heisenberg model. 

To construct the pure infinite-volume ground states, the symmetry breaking of two types 
discussed in the above should occur; A number of low-lying states of both types, a and j3, 
become degenerate to the ground state and rearrangements occur between them. This is con- 
sistent with the previous arguments by Bernu et a/.©, and Koma and TasakiQ. 

By discussing the symmetry breaking of two types separately, we find that the spontaneous 
sublattice magnetization is created by rearrangements of the low-lying states of the same type 
and that the spontaneous chirality is by rearrangements of pairs of low-lying states of the type 
a and of f3. 

It is possible that, if long-range order of one type exists in the true system, only one type of 
symmetry breaking occurs in the thermodynamic limit. Though the result from the renormal- 
ization group© indicates existence of one critical point (critical spin S c ), our arguments suggest 
that symmetry breaking can occur separately and that an intermediate phase can appear in 
which symmetry breaking of only one type occurs. 

An approximation for this state is given by 

1 4) = C/|F-order) (3.10) 
and the expectation values are related with the quantities of the ferromagnet in the forms 

(4\Uj2 S f u H 4 ) = (F-order| ^Sf |F-order) (3.11) 

i i 

and 

(4|Q Z |4) = (2\Q Z \2). (3.12) 



4 Evidence of symmetry breaking 

In this section we numerically study low-lying states of finite systems to verify the occurrence of 
symmetry breaking which was discussed in Section |3|. We display an evidence of the symmetry 
breaking in the Heisenberg and XY antiferromagnets. 

First we discuss the necessary conditions for symmetry breaking to occur. Bernu et al. 
have also discussed these matters. We discuss the conditions of the creation of the spontaneous 
sublattice magnetization and of the spontaneous chirality separately. 

When the symmetry breaking occurs in the thermodynamic limit, the low- lying states which 
satisfy the following two conditions should merge into the ground states^ in the infinite- volume 
limit. 



For a definition of the ground state in quantum spin systems, see refs. 23 and 1 
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(1) The energy (per site) is the same as that of the ground state in the infinite-size limit. 

(2) The spatial-symmetry is the same as that of the ground state of the classical model, i.e. 
the state has sublattice-translational invariance and Csy-invariance. 

In the thermodynamic limit, these low-lying states form pure ground states. In general, all the 
pure ground states give the same physical quantities. The sublattice magnetization and the 
chirality per site should be the same throughout the ground states. From this consideration 
and the discussions in Section 3, we obtain the necessary condition of the creation of the 
spontaneous sublattice magnetization as follows: 

(3) A number of low-lying states of the type a satisfy the condition (1) and they have the 
same macroscopic value of long-range order of the sublattice magnetization. 

The necessary condition of the creation of the spontaneous chirality is as follows: 

(4) Pairs of low-lying states of the type a and of (3 satisfy the condition (1) and they have 
the same macroscopic value of long-range order of the chirality. 

Bernu et al. studied the low-lying states of the Heisenberg model on the triangular lattice. 
They found many low-lying states which satisfy the conditions (1) and (2), and found that these 
states have similar macroscopic values for the sublattice magnetization. They estimated the 
sublattice magnetization of the ground state in the infinite-volume limit for the Heisenberg 
model. Leung and Runge® estimated the sublattice magnetization and chirality of the ground 
state for the XY and Heisenberg models. 

To examine these conditions explicitly, we study the energy gaps, the sublattice magneti- 
zation and the chirality of the low-lying states in the S = 1/2 XFand Heisenberg models, and 
estimate the physical values of the low-lying states in the infinite- volume limit. The systems 
of the size N = 9, 12, 21, and 27 with periodic-boundary conditions are studied, using the 
exact-diagonalization method. We study the low-lying states which belong to the subspaces of 
Sf otal = (or 0.5) and of 5f otal = 1 (or 1.5), and which satisfy the condition (2). (The ground 
state and some low- lying excited states really belong to these subspaces.) As we have indicated 
in Section 2, there exist the C 6 -symmetric (type a) and C 6 -antisymmetric (type (3) states. In 
many cases, we studied the lowest state of each subspace. There are exceptions: In the space 
of S'totai = for iV = 12 and Sf otal = 1.5 for N = 21 of the Heisenberg model, we chose the 
first-excited state of the type (3, since the first-excited state of the type a and the lowest state 
of the type (3 are degenerate and form a paired doublet, and the ground state (the type a) and 
the first excited state of the type f3 are in pairs. 

The lowest state in the S* otal = and C 6 -symmetric (type a) subspace is the ground state, 
which has been already calculated up to N = 36 in refs. [1^ and §. We used the data for iV = 36 
in ref. 

First we discuss the energy gap. Here the energy gap is defined as the difference between the 
total energies of two states. The energy gap between the lowest state of the type a and that of 
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(3 is shown in Fig. 2. In many cases the states of the type a and (3 are degenerate. The energy 
gap between the ground state (the type a) and the lowest state (the type a) of the Sf ota i = 1 
subspace is plotted in Fig. 3. It decreases proportionally to N~ x . As we have discussed in 
Section the former energy gap is relevant to the symmetry breaking of the chirality and the 
latter is to the breakdown of the rotational symmetry which creates the spontaneous sublattice 
magnetization. The latter energy gap corresponds to the singlet-triplet gap in the Heisenberg 
model, which was reported in ref. [5| 

It was proved that, if there is a long-range order, the energy gap decreases in the N~ l form 
or faster than it.©!^ Thus the above calculated results indicate an evidence for the existence 
of long-range orders in the XY and Heisenberg antiferromagnets. 

Next we consider the sublattice magnetization which is observed with the operator 

M = J2 S ? + E("^f + + £(4s? - ^Sf). (4.1) 

We calculated the long-range order of the sublattice magnetization (M 2 )/N 2 . The values of the 



spontaneous sublattice magnetization m are estimated, using the relation m = limjv-j-oo J6(M 2 ) /N 



for the Heisenberg model and m = limjvfoo 2^/ (M 2 )/N for the XY model. We show a deriva- 
tion of these relations in Appendix B. We fit the data of the sublattice magnetization m to the 
N^ 1 / 2 form. This finite-size correction was derived from the spin- wave theory© and using the 
effective Hamiltonian for a large spin©. (We also fit the long-range correlation (M 2 )/N 2 to 
the N~ l l 2 form, using the data for N = 3, 9, 12, 21, and 27, and found that the correction 
term of 0{N~ X ) in the fitting of (M 2 )/N 2 is larger than that in the fitting of m. Therefore we 
adopted the fitting of m.) The results for the XY model are shown in Fig. 4. We extrapolated 
in the N~ l l 2 form, using the data for iV = 9, 12, 21, and 27. For the ground state we also used 
the data for N = 36 which was reported in ref. |^. The values in the infinite-volume limit are 
shown in Table |[ They coincide with each other. Thus the condition (3) is satisfied in the XY 
model. The data for the Heisenberg model are shown in Fig. 5. The values in the infinite-size 
limit are shown in Table [| They are almost equal. Thus the condition (3) is also satisfied in 
the Heisenberg model. All the values are nonvanishing and close to the result by the spin-wave 
theory. However, as it is reported in ref. || if we fit the data (M 2 )/N 2 to iV -1 / 2 , we obtain 
the estimate m ~ 0. The values seriously depend on the fitting form. Thus it is still hard to 
conclude the existence of the long-range order in the Heisenberg model. 
Finally we consider the chiral order, which is observed with the operator 

SiXSj, (4.2) 

where the symbol % — ► j goes from the sublattice A to B, B to C, and C to A. We calculated 
the long-range order of the chirality ((Q z ) 2 )/N 2 . The values of the spontaneous chirality q 



are estimated using the relations q = lim^too y 3{(Q Z ) 2 } /N for the Heisenberg model and q = 
limAr-foo J ((Q z ) 2 )/N for the XY model. We show the results for the XFmodel in Fig. 6 and for 



12 



Table 2: Estimates of the sublattice magnetization m and the chirality q for the S = 1/2-XY 
antiferromagnet on the triangular lattice. 







m 


q 


^ 2 — n 

° total - u ' 


Type a 


0.41 


0.76 


°total — U ' 


Type (3 


0.39 


0.76 


Cz _1 
°total — 


Type a 


0.42 


0.77 


° total - x ) 


Type /3 


0.41 


0.78 



Spin-wave expansion 0.437® 0.798® 



the Heisenberg model in Fig. 7. For the XY model we extrapolated the expectation values of 
the chirality in the A r_3//2 form, which we derive from the finite-size correction of the spin-wave 
theory®. On the other hand, in the Heisenberg model the finite-size correction behaves in the 
iV -1 / 2 form, as discussed by Azaria et al. This difference of the correction terms comes from 
the fact that in the Heisenberg model the chirality is sensitive to spin-wave fluctuations of long 
wave-length, while in the XY model the chiral order is stable against them. The extrapolated 
values are shown in Tables |2| and |||. In the XY model the values are consistent with each other 
and nonvanishing, which suggests the existence of the chiral order. Thus the condition (4) is 
satisfied. In the Heisenberg model the values do not agree with each other. It is necessary to 
calculate larger systems to conclude that these low-lying states have the same chirality. We 
can not fit the data to a form including higher-order terms, since the size of the data is too 
small. Thus we can not still conclude the existence of the long-range order of the chirality in 
the Heisenberg model. 

To test the mechanism of symmetry-breaking of the chirality which we have discussed in 

Table 3: Estimates of the sublattice magnetization m and the chirality q for the S = 1/2- 
Heisenberg antiferromagnet on the triangular lattice. 





m 


q 


Stotai = °> T yP e a 


0.24 


0.21 


Stotai = 0, Type (3 


0.22 


0.16 


^totai = 1, T yP e « 


0.24 


0.54 


Stotai = X > T yP e P 


0.38 


0.73 


Spin-wave expansion 


0.250® 


0.405® 
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Section |3|, we construct the following state 

|0> = ^(|a)+i|/3», (4.3) 

where \a) (|/3)) denotes the lowest state of the type a (/3) of each S* otal subspace. We observe 
the spontaneous chirality of this state. The expectation value of the chirality is calculated 
through the relation 

(<j ) \Q z \<P)=i(a\Q z \P). (4.4) 

The estimates in the XY model are shown in Fig. 8. They are almost equal to the values 
which are shown in Fig. 6. This coincidence suggests correctness of the mechanism of the chiral 
symmetry breaking which was discussed in Section |3|. In the Heisenberg model we multiply the 
factor \f?> to the calculated values, since the symmetry is not fully broken in the state (|4.3| ); 
the state ( |3.7| ) should be properly used instead of (|4.3j). The estimates are shown in Fig. 9. 
Thus we can break the chiral symmetry by making a linear combination of low-lying states of 
the type a and of the type (3. 

From the above results we summarize that the necessary conditions (3) and (4) are satisfied 
in the XFmodel. In the Heisenberg model the necessary conditions are almost satisfied, though 
there still remains ambiguity in the point that the low-lying states have the same properties. In 
any case, in the above results there is no evidence which contradicts the occurrence of symmetry 
breaking. 



5 Summary 

To summarize, we studied the low-lying states of quantum antiferromagnets on a triangular 
lattice to clarify what kind of roles they play in the symmetry breaking and to show an evidence 
of occurrence of the symmetry breaking. 

We gave approximate forms of low-lying states, which have two-fold structures. It was found 
that these approximate states resemble the true states in various properties. We classified the 
low-lying states by the eigenvalues of C 6 , namely C 6 -symmetric and C 6 -antisymmetric. States 
of both types exist in pairs in the true low-lying spectrum. We discussed how rearrangements 
of the low-lying states of two types bring on the symmetry breaking. The spontaneous chirality 
is created by taking a linear combination of pairs of low-lying state of the type a and of j3. The 
spontaneous sublattice magnetization is obtained by making a linear combination of low-lying 
states of the same type. 

To display an evidence of the occurrence of symmetry breaking, we studied low-lying states 
of finite systems using the exact-diagonalization method. We studied whether the necessary 
conditions of the symmetry breaking are satisfied or not. We found that the conditions are 
satisfied in the XY antiferromagnet. The estimate of the sublattice magnetization is m = 0.41 
and of chirality is q = 0.77. These values are less than the spin- wave results by 7% and 4% 
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respectively. In the Heisenberg antiferromagnet the conditions are almost satisfied, although 
the estimates of the sublattice magnetization and the chirality do not converge well for lack 
of the system-size. The sublattice magnetization is estimated as m = 0.22 ~ 0.38 and the 
chirality is q = 0.16 ~ 0.73, where the results from the spin-wave expansion are m = 0.25 and 
q = 0.45. Since the extrapolated results seriously depend on fitting forms, it is still hard in 
the present study to conclude definitely the existence of the long-range order in the Heisenberg 
model. 



Appendix A. Low-lying states of the ferromagnetic XXZ 
model 

We show the uniqueness of the ground state and give spatial quantum numbers of the low-lying 
states of the spin S XXZ ferromagnet on finite systems. We discuss the model on the finite- 
volume triangular lattice A with a periodic boundary condition. For simplicity, we restrict the 
number of sites to an even integer. The Hamiltonian is given by 

h f = -j J2 ( s t s ! • si's 1 ; ■ //.s;.s;). (A.i) 

(ij>6A 

where J > and rj < 1. The summation runs over all the nearest-neighbor sites. 

The Hamiltonian (|A.1|) has nonpositive off-diagonal elements. The basis which have the 
same S^ otal are connected by the elements of the Hamiltonian. From the Perron- Frobenius the- 
orem, the lowest eigenstate in each S* otaX subspace is unique and it has nonnegative coefficients. 

Here we show that the ground state exists uniquely in the S* otgl = subspace. We consider 
the following Hamiltonian 

H' F = ~ E MSfSf + S?S] + r}S*S*), (A.2) 
ijeA 

where Jjj > for all i and j. Affleck and Lieb© showed that, in the antiferro magnetic XXZ 
chain, the ground state exists uniquely in the S* otal = subspace. As they did, we transform 
the Hamiltonian ( |A.2j) using the unitary operator 



H/ = exp(^£ 5 f) (A.3) 

into the form 

W^H' F W = - J £ {1(1 + n)(S+S- + SrS+) + 1(1 - V )(S?Sf + SrSr) + S*S*}. (A.4) 

i,jeA 

In the region — 1 < rj < 1, this transformed Hamiltonian has nonpositive off-diagonal elements. 
It has two connected blocks, i.e. all basis in the subspaces with even 5f otal are connected by 
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the elements of the Hamiltonian and those with odd Sf tai are connected. From the Perron- 
Frobenius theorem, the lowest eigenstate in each connected subspace is unique. Thus we find 
that the number of the ground states is at most two in the original model (|A.2|) . The lowest 
states in the S* otal = n and Sf otal = — n subspaces are degenerate. When level crossing of 
the ground state occurs, more than two ground states must exist, which contradicts the above 
result. Thus in the whole region of various parameters { Jij }, the ground state should exist in 
the same S* otgl subspace. 

For the region 77 < —1, as Affleck and Lieb did, the above results can be shown using the 
unitary operator 

W = exp(£ y £S?). (A.5) 

To obtain the eigenvalue S* otal of the ground state, we consider the case = J for all i 
and j. The Hamiltonian can be written in the form 

H' F = -J{(Sy-(l- V )(Srh (A.6) 

where S= (S x ,S y ,S z ) and S a = ^2i e \S" (a = x,y,z). The ground state of this model is 
unique and it has the eigenvalues S to tai = SN and Sf otal = 0.® From the above results, we 
find that the ground state of the original model ( |A. 1|) in the regions 77 < — 1 and — 1 < 77 < 1 
exists uniquely in the S^ otal = subspace. 

This argument can be extended to the ferromagnetic XX Z model on any lattice. The ferro- 
magnet on a bipartite lattice is equivalent to the XX Z antiferromagnet with the ^-component 
anisotropy —rj. In this case the above results are consistent with what were proved by Lieb and 
Mattis©, and Affleck and Lieb©. 

Lastly, we give the spatial quantum numbers of the lowest state in each S* otal subspace. 
The Hamiltonian ( |A.1|) is invariant under any translation, rotation, and reflection. As we have 
shown in the above, the lowest state is nondegenerate and it has nonnegative coefficients. From 
this fact, we find that the lowest state has the eigenvalue 1 for any translation, rotation, and 
reflection, i.e. it is translationally invariant and Cgy-invariant. 

Appendix B. Relations between long-range order and 
symmetry breaking 

We discuss the relation between spontaneous symmetry breaking m and the long-range order 
parameter a. It has been discussed® 0) that the relation is given by m = \f?>a for the 
Heisenberg antiferromagnet and by m = \p2o for the XY antiferromagnet on bipartite lattices. 
But these relations are not valid in the antiferromagnets on the triangular lattice. Here we 
show that, in the antiferromagnets on the triangular lattice, the factor becomes a/2 times as 
large as the above. 
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Our arguments are based on the assumption that the symmetric infinite-volume ground 
state, ( |3.1| ), is decomposed into the pure ground states. Koma and Tasaki® used this decom- 
position to explain the factor v3 for the Heisenberg antiferromagnets on bipartite systems. 

The spontaneous sublattice magnetization m is defined as 

,. ,. lTr[Mexp{-/3(H-BM)}] 
m = hm hm hm [ - , l , v B.l) 

bio tvtoo /3too N Tr[exp{—(3(H — BM)}] V ; 

Using the sublattice-translational invariance and equivalence of the sublattice magnetization, 
we have 

m = uiS*), (B.2) 

where the state uj{- ■ ■) is defined by ( |3.5| ). The long-range order parameter of the sublattice 
magnetization a m is defined by 

tvtoo moo N 2 Tr[exp(-/3#)] v ; 

Using the rotational and translational invariance, we obtain 

^ 2 =limi { u,(S* S?)+u(-~S* S? +ei + y^S x S v r+ei ) 



+u(-oSSS^ - V^+e 2 )}, (B.4) 



1 \/S 

_ ox ox v . 

2°0 O r+e 2 2 

where the state uj is defined by Q3.1] ). The sites and r belong to the sublattice A and ei (e 2 ) 
denotes the unit lattice vector, which belongs to the sublattice B (C). 
We also define the spontaneous chirality q as 

, , , 1 Ti\Q z exv{-(3(H - BM)}] NN , , 

g = limlim lim— — ^ f 1 n ,) T — — ^ = u(Q z (0)) B.5 

H bio tvtoo pioo N Tr [exp{- [3 (H - BM)}} ^ i v > 

and the long-range order parameter of the chirality a q as 

a q 2 = hm hm ± T^l^L^)l = l im w (Q*( )Q*( r )), (B.6) 

9 TVToo /3Too AT2 Tr[eX P (-/3iy)] |r|Too ^ v y ^ v v ; 

where the operator Q 2 (V) is the z-component of the chiral order-parameter operator on the unit 
triangular cell, 

2 

From now we derive the relations between m and <r m , and between q and cr g . Our results 
are as follows: For the XY model 

m = 2a m and q = <J q , (B.8) 

and for the Heisenberg model 

m = VQcr m and q = V3a q . (B.9) 
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The relation m = \/Ga m for the Heisenberg model has been already used by Bernu et a/.©. 

To derive these relations, we use the following standard arguments. It is widely believed 
that the mixed state u can be naturally decomposed into pure equilibrium states,© 

oj{...) = f dau a {---), (B.10) 



where {u a } denote the pure states and the parameter a describes the properties of them. The 
state uj, which is defined in eq. ( |3.5| ), is one of u a . The pure states have the cluster property,® 

uj a (A B r ) — >u a (A )u a (B r ) (B.ll) 

|r|foo 

for any local operators A and B. 
B.l XY-like model 

In the XXZ model with —0.5 < A < 1, it is expected that the decomposition (|B.10|) is of the 



form 

UJ 



) = Jl f^dOicu ,,+(.-.) + w 9,-{"-)h (B.12) 

47T JO 

where ojq i+ (uq-) denote the states which have positive (negative) chirality and in which vectors 
of spins on the sublattice A form the angle 9 with the x-axis. The state u corresponds to 
Ug=Q +. The pure states have following expectation values of the sublattice magnetization 

u$ j+ (sfj — tocos 6*, ue,-\Sf) — tocos 0, (B.13) 

^, + (-\s x ei + ^S* ) = mcos9, uoe,-(-\s x ei + ^-Sl x ) = mcos(9 - 

"e, + (-\s: 2 - ^Sl 2 ) = mcostf, ue,-(-\si 2 - ) = mcos(9 + 

and the chirality 

uj e ,±(Q Z (0)) = ±q. (B.14) 

Using the decomposition ( B.12|) and the property ([B.llj) , the long range order parameter 
of the sublattice magnetization, (|B.4| ), is transformed as 



a, 



1 f 2n r 1 \/3 

~~ iHP To - / ^{^-(^o^r) + V0,+(--S%S* + + —S$S? +fll ] 

|r|Too 1Z7T JO Z Z 



1 

1 /" 27r 



127T JO 
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= ^! f 2n d e C os 2 e = ^. (B.i5) 

An Jo 4: 
Thus we obtain the relation m = 2a m . 

In the same way, the long-range order parameter of the chirality, (|B.6| ), is estimated as 

IHtoo An Jo 



1 r 27T r 

a q 2 = lim-/ deU + {Q Z m z {r))+uJeAQ Z m z {r))} 

r Too 47T JO 1 J 



= -L r^H, + (g 2 (o)) 2 + ^,_(g 2 (o)) 2 } = g 2 . (B.ie) 

47T JO 

Thus we get the relation q = a q . 
B.2 Heisenberg model 

For the Heisenberg model it is expected that the mixed state uj{- ■ ■) is decomposed as 

">(•••) = A / r dc^uaA---), (B.17) 



57T 2 J JO 



where the vector f2 is perpendicular to the plane to which all spins are parallel, and the angle 
4> denotes rotation with respect to the vector Q. The direction of the vector f2 and that of 
the chirality are the same, which are described with the spherical coordinates (9, (p). The state 
Co corresponds to the state u v =o. The pure states have following expectation values of the 
sublattice magnetization 



u!n,4> ySy) = m(cos (f cos 9 cos — sin cp sin <p) , (B. 18) 



2 ei 2 

mjcos^o — — ) cos cos(^0 + — ) — sin^y? — — ) sin(0 + — ) |, 

V5. 



= mjcos^o + — ) cos cos (0 — — ) — sva(ip + — ) sin^0 — §")}' 
and the chirality 

w n ,*(Q(0)) = gO. (B.19) 
Using the decomposition ( B.17|) and the property (|B.11|) , the long-range order parameter 



of the sublattice magnetization is transformed as 



2 - l^Jdn fj d^n^Sl 



2 

m 



6 (R20) 
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Thus we obtain the relation m = \^Qa m . 

Using the decomposition ( |B.17| ) and the property ( B.ll ), the long-range order parameter 
of the chirality is transformed as 

a q 2 = i lim u(Q(0)Q(r)) (B.21) 

O \r\Joo 

q 2 r r2ir 



24vr 2 

where the isotropy of the state is used. Thus we get the relation q = V3cr q . 
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Figure Captions 



Figure 1: Two-fold structures in the spectrum of low-lying states. 

Figure 2: Size dependence of the energy gap between the lowest state (the type a) and the 
lowest excited state (the type (5) in the same S' t z otal -subspace for the Heisenberg model. 
Figure 3: Size dependence of the energy gap between the ground state (the type a) and the 
lowest state (the type a) in S^ otal =l. 

Figure 4: Size dependence of the sublattice magnetization of the lowest state in each subspace 
of the XY model, which is estimated through 2^J (M 2 ) /N. The data for iV=36 is also listed©. 
SW denotes the result from the spin- wave expansion®. 

Figure 5: Size dependence of the sublattice magnetization of the lowest state in each subspace 
of the Heisenberg model, which is estimated through yj6(M 2 ) /N. The data for N=36 is also 
listed©. SW denotes the result from the spin- wave expansion w. 

Figure 6: Size dependence of the chirality of the lowest state in each subspace of the XY 
model which is estimated through yj ((Q z ) 2 ) / N . The data for iV = 36 is also listed®. SW 
denotes the result from the spin- wave expansion®. 

Figure 7: Size dependence of the chirality of the lowest state in each subspace of the Heisenberg 
model which is estimated through ^3((Q Z ) 2 )/N. The data for iV=36 is also listed®. SW 
denotes the result from the spin- wave expansion®. 

Figure 8: Spontaneous symmetry breaking of the chirality of the XY model which is created 
by a linear combination ^(|a) + of the low-lying state of the type a and that of (3. 
Figure 9: Spontaneous symmetry breaking of the chirality of the Heisenberg model which is 
created by a linear combination -^(\ce) + of the low-lying state of the type a and that of 
[3. The values are multiplied by the factor \^3. 
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